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Abstract 

A sufficient condition for tfie linear stability of three dimensional equilibria 
with incompressible flows parallel to the magnetic field is derived. The con- 
dition involves physically interpretable terms related to the magnetic shear 
and the flow shear. 
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I. Introduction 



For static ideal magnetohydrodynamic (MHD) equilibria there is a pow- 
erful tool known as "the energy principle" providing necessary and sufficient 
conditions for linear stability [1]. In the presence of flow, however, the sta- 
bility problem is much tougher because the force operator becomes non Her- 
mitian; thus, only sufficient conditions were obtained [2]-[TU]. Motivation of 
the present study is a couple of papers by Ilin and Vladimirov [SI [S] in which 
a sufficient condition was derived for the linear stability of plasmas with con- 
stant density and incompressible flows parallel to the magnetic field. This 
condition states that an equilibrium is stable to three dimensional pertur- 
bations provided that: i) the flow is sub-Alfvenic and ii) inequalities (51) of 
Ref. [H] are satisfied. Here we show, however, that those inequalities are not 
correct for the following reasons: 

1. The authors of Refs. [B] and [S] have not noticed that because of 
the field aligned flow the equilibrium current density lies on magnetic 
surfaces. This property simplifies the stability analysis and results in 
a single inequality for the sufficient condition in place of the couple of 
inequalities (51) of Ref. [S]. 

2. A term associated with the flow shear was ignored in Refs [6] and [8]. 

The correct sufficient condition obtained here contains physically interpretable 
terms related to the magnetic shear and the flow shear. 

The equilibrium characteristics are examined in Sec. II including a prove 
of the coincidence of the current density surfaces with the magnetic surfaces. 
Sec. Ill reviews the energy principle established in Refs. [B] and [5] which 
subsequently is employed in Sec. IV to derive the sufficient condition. A 
major part of the derivation is presented in the Appendix. 

II. Equilibrium 

We consider the steady states of a plasma of constant density and incom- 
pressible flow parallel to the magnetic ffied in the framework of ideal MHD 
(see for example Eqs. (l)-(6) of Ref. [11] written in convenient units and the 
density set to unity). Also, it is assumed the existence of well defined equi- 
librium magnetic surfaces in three dimensional geometry which are labeled 
by a smooth function ip. Using 

V = AB, (1) 
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where A is an arbitrary function, the incompressibihty condition (V ■ V = 0) 
imphes that A is a surface quantity: 



A = A(V^). (2) 

Then, employing the identity (V ■ V) V = Vy^/2 - V X V X V, the momen- 
tum equation 

(V- V)V = JxB- VP 

leads to 

(l - A^) JXB = V (P + ^) - P^V , (3) 

where P is the magnetic field modulus. The component of along the 
magnetic field implies that the quantity P + A^P^/2 is uniform on magnetic 
surfaces: 

A^P^ 

P+^^PsW. (4) 

Thus, owing to the flow the isobaric surfaces depart from the magnetic sur- 
faces unlike the case of static equilibrium associated with the surface function 
Ps{ip)- Consequently, Eq. ([3]) is put in the form 

1 - A^) JxB = Pyij - (A^)'— 
or 

N = JxB = ^(V^,p2)Vz/; (5) 

where 

Eq. (ISl) implies that the current density lies on magnetic surfaces a prop- 
erty not noticed in Refs. [6] and [8]. Note that this holds because of the 
incompressible field aligned flows; for flows of arbitrary direction the current 
surfaces do not coincide with the magnetic surfaces. The fact that B, J and 
V share the same surfaces simplify the stability analysis to follow. To this 
end we also will need the quantity 

M = VxN = V^xV^ (7) 

from which it follows that 

M ■ N = 0. (8) 
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III. Review of the energy principle 

In Refs. [HI [H] an energy principle was established for incompressible 
perturbations [V ■ ^(x, t)=0] around a steady state and non-slip boundary 
conditions: 

V ■ n = b ■ n = ^ ■ n = 0. (9) 

Here v(x, t) and b(x, t) are the perturbations of the velocity and the magnetic 
field and conditions are imposed on a fixed boundary &D surrounding the 
plasma domain V. The principle is based on the fact that the perturbation 
energy 

[ f if - 1^ ■ dV, (10) 



V V2' 2 

is conserved by the linearized ideal MHD equations {dE/dt = 0). Here K is 
a symmetric operator defined by the formula 

A'^ = VxVxv + vxi7-BxVxb-bxJ, 

where v = Vx (^xV), 

b = Vx(|xB), (11) 

and ft = Vx V. Evidently, E a.s a. quadratic functional of ^ and ^ is positive 
definite if the potential energy 

w = -\l i- k^dv (12) 

z Jv 

is positive definite. It is known, however, that for flows of arbitrary direction 
the functional W is never strictly positive definite [2]-[T0]. For this reason 
further consideration is restricted to the steady states with field aligned flows 
described in Sec. II. In this case ( !T2|) can be written in the form 

1^ = 1 {(1 - A^) [b^ + b ■ (JxO] - 2A(^ ■ VA) ■ (B • V)B]]dV. (13) 

Derivation of f[T^ is given in Ref. [B]. Whenever the potential energy (1131) 
is positive definite the equilibrium is linearly stable. 

IV. Sufficient condition for linear stability 
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As in Refs. [H] and [S] assuming that J X B 7^ we express the perturba- 
tion vector ^ in the form 

C = a(x,t)N + /5(x,t)J + 7(x,t)B. (14) 

It can then be shown (see Appendix) that W assumes the form 

W = Wi + W2, (15) 

Wi = l- f {1- X'^){h + a3xNfdV, (16) 
2 Jv 

W2= f Aa\ (17) 



where 



/l = -(l-A^)(JxN)-(B- V)N-A(N-VA) IN-^— + ivM. (18) 

Evidently, W is positive semidefinite if |A| < 1 and 

A>0 in V. (19) 

Inequahty (|T9l) is substantially different from the respective inequalities (51) 
of Ref. [8]. In particular, the last term of (fT8|) containing VA was missed in 
[6] and [8]. Using the equilibrium relations of Sec. Ill, (fT8|) reduces to 

A = -g'^{{l-Xf{3xVtlj)-{B-V)Vtp 

-^-^mf(v^.^ + g\V^A]. (20) 



On account of (|T5 i) -( fT7l) and ( l20l) we can conclude that a general steady 
state of a plasma of constant density and incompressible flows parallel to the 
magnetic is stable to small three-dimensional perturbations if i) the flow is 
sub-Alfvenic and ii) 

~. A 

A = — >Q. 21 

9^ 

Using the relation 

(B ■ V)V^ = JxV^ - (VtA ■ V)B 
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A can be put in the physically interpretable form: 



A 



(1 - A)2 {JxV^Y - (Jx V^) ■ (VV^ ■ V)B 



(22) 
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The first negative destabilizing term in (122!) should be related to current 
driven modes. The other terms can be either stabilizing or destabilizing. 
This depends on the sign of (A^)' in relation to the velocity shear and on the 
differential variation of B and perpendicular to the magnetic surfaces in 
relation to the magnetic shear. Also, the last term has an additional implicit 
dependence on (A^)' and P'^ through the quantity g [Eq. ([6])]. 

It is recalled that the sufficient condition established here can be applied 
to any steady state without geometrical restriction. Application to steady 
states of fusion concern in connection with possible stabilizing effects of the 
flow is under way. 



The procedure to follow is based on that of Appendix of Ref . [8] . Since 
there are substantial differences, however, the derivation will be presented in 
a self contained way. 

Preliminarily, in view of the representation f|T^ for ^ and the incompress- 
ibility condition V ■ ^ = we obtain the following relations: 



Appendix: Derivation of (ll5l)-(fT8l) 



V ■ (aN) + J ■ V/3 + B ■ V7 = 0, 



(23) 



^xB = aNxB + /?N, Jx| = aJxN + 7N. 



Also, ffTTl) becomes 



b = Vx(aNxB + 6N) 



The first term of f|T3|) is written as 




(24) 
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Employing the identity V-(AxB) =A-VxB — B-VxAwe have in 
connection with the last term of 



b-(7N) = 7N- Vx (aNxB + /?N) 

= (aNxB + /3N)- Vx(7N) + V - [(aNxB + /3N) x (7N)] 



and therefore 

(1 - A2)b • (7N) 



(1 - A^) (aNxB + pN) ■ Vx(7N) 
+V ■ [(1 - A2)7(aNxB + /3N) xN 



the last term following from the fact that N is proportional to 'Vip [Eq. ([5]) 
Substituting fl25l) into flMl) and integrating by parts furnishes 



(l-A^) (b + aJxN)2 -a2(jxN) 



2Jv 



Y2 



-a(JxN) • b + Vx(7N) ■ (aNxB + /5N) 



dV. 



(25) 



Furthermore, employing (JTj) and (IHl) we find for the above defined quantities 
Yi and Y2: 

= Vx(7N) ■ (aNxB + /5N) = (7VXN + V7XN) ■ (aNxB + /3N) 
= a7M ■ NxB - aiV2(B ■ V7) 
Y2 = -a(JxN) ■ Vx (aNxB + /5N) 

= -a(JxN) ■ [Vax(NxB) + aVx(NxB) + V/3XN + /3M] 

= -a2(JxN) ■ Vx(NxB) -aAr2(N- Va + J- V/3) -a/3M- (JxN) 

and with the aid of (l23l) 



Yi+Y2 



-a 



(JxN) ■ Vx(NxB) - N^V ■ N -a/3J-(NxM)-a7B-(NxM) 



Eq. (125!) then becomes 

I = \ l (l-A2)f(b + aJxN)2 

(JxN) ■ Vx(NxB) - A^^V ■ N + (JxN)^ 



—a 



-a/3J ■ (NxM) - a7B ■ (NxM)} c/y. 



(26) 
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The staff within the square brackets in fl26p can be put in the concise form 
(JxN) ■ Vx(NxB) - A^^v . N + (JxN)^ = 2(JxN) ■ (B ■ V)N. (27) 
To show this we employ the relation V(B ■ N) = implying that 
(N- V)B = -(B- V)N + MxB + JxN; 

then, 

(JxN) ■ Vx(NxB) = (JxN) ■ [(B- V)N-B(V-N) - (N- V)B] 

= 2(JxN)-(B- V)N-(JxN)2 + iV2V-N. 

We now consider the second part of W [see Eq. f|T3|) ]: 

Q==- [ A(^ ■ VA) ■ (B ■ V)B] dV. (28) 
Using the relations 

(C ■ A) = aN ■ VA 

and 

^■(B- V)B = ^- (^- + JxBj = {aN + 133 + -fB) ■ + aN^, 



281) is put in the form 

Q = -/Ja„=(N.VA)(n.:^ 



N' 
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+Aa/3(N- VA)(J-^— ) + Aa7(N-VA) IB-^— J \ dV (29) 

The first term in fl29l) containing VA was ignored in Refs. [6] and [8]. In 
view of (ESI), (EZD and (1291), is written in the form 



W = Wi + W2, (30) 

Wi = -f{l-X^)(h + a3xNfdV, (31) 
2 Jv 
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Wo 



1 

2Jv 



[-2a^ [(1 - A2)(JxN) ■ (B • V)N 
/ 



+A(N ■ VA) N 



—al3 



-07 



;i - A2)J ■ (NxM) + 2A(N ■ VA) (^J ■ 
(1 - A2)B ■ (NxM) + 2A(N ■ VA) ( B 



Vfi2 



dV. 



(32) 



The coefficients of a(3 and 07 vanish identically. Indeed, on account of the 
equilibrium relations ([2]), (jlj), ([5]) and ([6]) we have for the coefficient of a(3: 



[I - A2)J ■ (NxM) + 2A(N ■ VA) J 



Vfi2 



:i-A^)(?|V^|^J-V(7+(A^)'^7|V^|^J 



V52 



J V 



R2 



0. 



The coefficient of 07 also vanishes because is symmetric to the coefficient of 
af3 in replacing J with B. Consequently ( l32l) assumes the form ( JT71) . 
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